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SECTION  1 
INTRODUCTION 


Ionospheric  disturbances  have  been  observed  following  large, 
low-altitude  nuclear  weapons  tests  conducted  by  the  U.S.A.  in  1962l_2f 
the  Soviet  Union  in  1961 3  and  F ranee  in  19674.  Close  inspection  of  the 
data  reveals  that  there  appear  to  be  two  types  of  waves:  a  wave  which 
propagates  to  at  least  several  thousand  kilometers  from  the  hurst  and 
whose  period  depends  on  the  observation  point;  a  wave  which  propagates  all 
the  way  around  the  earth  (possibly  more  than  once)  and  whose  period  is 
around  15  minutes.  There  was  agreement  very  soon  after  the  observations 
were  made  that  the  Ionospheric  disturbances  were  due  to  a  system  of 
acoustic  gravity  waves  set  up  by  the  explosion. 

The  ionospheric  observations  inspired  a  substantial  research 
program  on  the  subject  of  acoustic  gravity  waves5.  During  this  effort 
much  was  learned  about  the  rather  unusual  features  of  the  propagation  of 
such  waves.  Of  particular  importance  was  the  classification6  of  the  waves 
as  either  freely  propagating  (which  can  exist  in  any  stratified  atmos¬ 
phere)  or  ducted  (which  can  exist  only  In  an  atmosphere  with  a  certain 
type  of  temperature  structure);  the  earth's  atmosphere  Is  such  as  to 
support  both  freely  propagating  and  ducted  types  of  waves. 

In  spite  of  the  progress  made  In  understanding  acoustic  gravity 
waves  the  fundamental  Issues  associated  with  the  test  data  were  never 
resolved.  In  particular  It  was  never  determined  what  aspect  of  the  explo¬ 
sion  Is  responsible  for  creating  the  gravity  waves  or  what  type  of  wave 
(ducted  or  freely  propagating)  was  most  Important  for  producing  the  iono¬ 
spheric  disturbances.  The  most  commonly  suggested  source  mechanisms  were 
the  blast  wave  and  motions  associated  with  the  rising  fireball.3’9  More 


recently  we  have  suggested  that  the  heated  region  above  100  km  altitude 
which  is  created  directly  over  a  low  altitude  explosion  by  the  restrength¬ 
ening  of  the  blast  wave  might  act  as  a  source  of  acoustic  gravity  waves.7 

In  1982  DARPA/AFOSR  tasked  MRC  to  use  mathematical  procedures 
developed  at  MRC  during  the  1970's  to  perform  calculations  to  determine 
what  the  most  Important  source  or  sources  for  the  nuclear-explosion-induc¬ 
ed  gravity  waves  are.  This  document  is  the  final  report  on  that  effort. 

We  have  made  calculations  for  the  coupling  of  the  blast  wave,  the  rising 
fireball  and  the  high  altitude  hot  spot  to  both  ducted  and  freely  propa¬ 
gating  acoustic  gravity  waves.  Although  the  simplified  nature  of  our 
atmospheric  models  precludes  detailed  comparisons  between  the  calculations 
and  the  data,  we  can  compare  the  general  features  of  the  two.  The  results 
allow  us  to  draw  these  conclusions:  1)  the  major  source  of  the  large 
amplitude  Ionospheric  gravity  wave  observed  following  the  U.S.  and  Russian 
tests  (the  wave  with  periods  ranging  from  30  minutes  to  more  than  two 
hours)  is  the  blast  wave  and  the  energy  is  transmitted  by  the  freely 
propagating  modes;  2)  the  source  of  the  "15  minute  period  waves  that 
were  observed  to  go  all  the  way  around  the  world  is  the  high  altitude  hot 
spot  and  the  energy  Is  traveling  in  the  ducted  modes;  3)  the  motions 
associated  with  the  rising  fireball  give  a  smaller  signal  in  the  freely 
propagating  modes  than  the  blast  wave  and  a  smaller  signal  in  the  ducted 
modes  than  the  high  altitude  hot  spot;  4)  it  appears  that  ionospheric 
observations  will  be  of  limited  use  In  detecting  small  (<10kT)  low  alti¬ 
tude  detonations  at  large  distances  (>1000km)  from  the  burst. 

In  the  next  section  of  the  present  report  we  briefly  review  the 
mathematical  procedures  developed  In  References  6  and  7;  we  then  formu¬ 
late,  within  the  context  of  these  procedures,  the  problem  of  the  coupling 
of  the  blast  wave  to  the  freely  propagating  acoustic  gravity  waves  in  an 
Isothermal  atmosphere.  The  solution  Is  given  as  a  sum  of  Integrals  which 
represent.  In  our  way  of  thinking  about  the  problem,  a  sum  over 
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eigenstates.  We  have  evaluated  these  Integrals  numerically  and  the 
results  are  presented.  The  fact  that  our  current  calculations  are  for  an 
isothermal  atmosphere  whereas  the  data  are  for  a  region  of  the  real  atmos¬ 
phere  over  which  the  temperature  varies  considerably  causes  certain  ambi¬ 
guities  In  comparing  our  results  with  data;  the  way  we  have  resolved  these 
ambiguities  is  described  in  detail  in  Section  2C.  The  results  of  the 
comparison  of  the  calculations  with  the  data  shows  good  agreement  for  both 
the  amplitude  of  the  disturbance  and  its  wave  form;  this  is  the  first 
calculation  of  the  acoustic  gravity  wave  generated  by  a  nuclear  explosion 
to  show  such  good  agreement.  We  would  expect  the  agreement  between  the 
calculations  and  the  data  to  be  improved  by  the  use  of  more  realistic 
atmospheric  models  (such  calculations  would  also  substantially  reduce  the 
ambiguities  mentioned  above)  but  the  agreement  is  already  sufficiently 
good  that  we  are  confident  that  the  blast  wave  is  the  source  mechanism  for 
the  large  amplitude,  varying  frequency  gravity  waves. 

Prior  to  resorting  to  numerical  procedures  for  evaluation  the 
integrals  we  attempted  to  estimate  them  in  the  stationary  phase  approxima¬ 
tion.  This  procedure  worked  well  enough  for  sufficiently  late  time  but 
the  stationary  phase  approximation  proved  surprisingly  ineffective  for 
times  during  the  first  several  oscillations  after  the  arrival  of  the 
signal  (the  interesting  times  for  our  work).  In  Section  2  and  Appendix  A 
we  give  the  details  of  our  stationary  phase  calculations  and  the  difficul¬ 
ties  encountered.  We  also  remark  on  the  relevance  of  the  results  to  the 
previous  work  of  other  workers.4 

In  Section  3  we  consider  the  coupling  of  the  motions  associated 
with  the  rising  fireball  to  acoustic  gravity  waves.  The  calculations  are 
done  with  relatively  simple  models  consisting  of  point  sources;  the  calcu¬ 
lations  are  done  in  the  approximation  that  only  low  frequencies  contri¬ 
bute.  The  reason  that  these  calculations  have  been  done  in  less  detail 
than  those  for  the  blast  wave  is  that  the  calculations  with  the  simple 


model  make  it  clear  that  the  motions  associated  with  the  rising  fireball 
are  a  much  weaker  source  of  gravity  waves  than  the  blast  wave. 

In  Reference  7  we  presented  calculations  which  gave  the  coupling 
of  the  high  altitude  hot  spot  to  ducted  gravity  modes.  In  that  report 
we  gave  specific  calculations  only  for  small  events;  the  formulas  given 
there  can  be  applied  to  larger  events  so  as  to  allow  comparison  with  the 
data  in  Reference  4.  The  comparison  between  the  calculation  and  the  data 
shows  good  agreement  for  the  amplitude  but  somewhat  less  good  agreement 
for  the  wave  form.  We  believe  the  differences  between  the  calculated 
waveform  and  the  observed  one  are  due  to  the  effects  of  dispersion  which 
is  absent  from  our  calculations  but  present  to  some  extent  in  the  real 
atmosphere.  Although  we  expect  the  dispersion  to  be  small,  for  the  ten 
thousand  kilometer  path  lengths  of  the  observations  the  waveform  could 
have  become  appreciably  distorted. 

By  combining  the  calculations  presented  in  this  report  with 
those  presented  in  Reference  7  we  arrive  at  the  conclusions  stated  above. 
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SECTION  2 

GRAVITY  WAVES  GENERATED  BY  THE  BLAST  WAVE* 


A.  INTRODUCTION 

In  this  section  we  investigate  the  gravity  wave  response  of  the 
atmosphere  to  a  low  altitude  nuclear  burst  using  methods  developed  earlier 
in  References  6  and  7.  Our  goal  in  this  section  is  to  determine  the  grav¬ 
ity  wave  signal  generated  by  the  blast  wave  and  carried  by  the  freely 
propagating  modes  and  to  compare  with  limited  experimental  observations1-5 
This  investigation  is  complicated  by  the  fact  that  the  Euler-Lagrange 
equations  for  a  stratified  atmosphere  are  nonlinear.  These  nonlinearities 
cause  lengthening  and  strengthening  of  the  initial  blast  wave  during  the 
early  time  development  (to  about  100  km  from  the  point  of  blast)  which 
continues  until  the  relative  overpressure  is  about  5  percent.  This  early 
time  development  has  been  described  by  a  model  in  Reference  8.  For  large 
distances,  say  1000  km  or  more,  the  relative  pressure  variation  is 
sufficiently  small  and  we  can  assume  that  the  hydrodynamics  is  well 
described  by  linearized  Euler-Lagrange  equations  in  an  isothermal  atmos¬ 
phere.  To  obtain  a  smooth  transition  from  the  nonlinear  regime  to  the 
linearized  one,  we  use  the  calculations  of  amplitude  A  and  pulse  length  L 
of  Reference  8  to  define  the  initial  boundary  condition  at  some  chosen 
height  D  for  the  linearized  differential  equations  for  the  atmosphere. 

From  this  point  on  the  linearized  equations  are  solved  numerically  using 
methods  of  References  6  and  7  and  these  results  are  compared  to  observa¬ 
tions  and  the  stationary  phase  approximation,  an  approximation  commonly 
used  in  estimating  gravity  wave  response.  As  we  shall  see  the  exact 
calculation  compare  favorably  with  available  observations  and  in  passing 

*  The  methods  used  in  this  report  were  developed  and  presented  in  Refer¬ 
ence  7.  For  completeness,  we  review  the  methods  in  this  Section;  much  of 
this  material  is  taken  from  Reference  7. 
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we  shall  note  that  the  stationary  approximation  is  very  inaccurate  in 
estimating  the  amplitude  of  the  gravity  waves  arising  from  a  blast  wave. 


B.  HYDRODYNAMIC  VARIABLES  AND  EQUATIONS 


To  begin  with,  we  consider  the  case  of  one  isothermal 
The  hydrostatic  solution  about  which  we  wish  to  perturb  is  given 


u  *  0 

P  *  P  e 
o 


-Z/H 


p  =  pQe 


-z/H 


layer. 

by 


where  u  is  the  fluid  velocity,  P  the  pressure,  and  p  the  density.  P0, 
p0,  and  H  are  constant  related  by  P0  =  PoHg,  where  g  is  the  acceler¬ 
ation  of  gravity.  We  define  perturbing  variable  by  the  relations 

S  »  V  ez'2H 
P  =  e-z/H(P0+P  ez/2H) 

p  -  e"z/H(p  +p  ez/2H). 
o 

We  have  chosen  the  exponential  factors  so  that  the  perturbing  functions 
V,  P,  p  will  be  bounded  for  all  values  of  x,y,z,  and  t,  and,  for  compact 
disturbances,  will  be  square  integrable  functions  of  x,y,  and  z  for  any 
value  of  t.  (This  fact  is  not  obvious  but  is  seen  In  the  results  below.) 
For  these  variables  the  linearized  hydrodynamic  equations  take  on  the 
familiar  form 


We  now  wish  to  construct  a  Hilbert  space  in  which  our  solution  will  lie. 
We  arrange  the  functions  into  a  column  vector  in  the  order 


and  define  the  inner  product 


00 

<a|b>  =  /  Pa  Pb  +  ^za  ^zb+  ^xa  ^xb  +  ^ya  ^yb  +  papb  dxdydz  (2-3) 


Equation  2  takes  on  the  form 

M j a >  3  <i)|a>,  (2-4) 

where  the  operator  M  is  given  by 
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(2-7) 


(2-8) 
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In  the  new  space  our  eigenvalue  problem.  Equation  (2-2),  becomes 


M I a>  =  <*la>  (2-9) 


where 


This  operator  is  indeed  selfadjoint  and  we  can  use  all  the  standard 
procedures.  If,  for  example,  the  perturbed  configuration  at  time  0  is  |f>. 
At  later  times  the  configuration  will  be* 

|a>  a  /<(u|f>  j tu>  e  ^^dw  (2-11) 

where  the  U>  are  the  eigenvectors  of  M  with  eigenvalue  w. 


t  Equation  2-11  symbolically  Indicates  that  we  must  sum  over  all  the 
eigenstates;  due  to  degeneracy  the  sum  will  be  a  multiple  Integral. 

10 


c. 


EIGENVECTORS 


Direct  substitution  into  (2-10)  readily  verifies  that7 
I 

Pi  Jo(kr) 


k,a>.  = 


p2  Jo(kr) 
P  3  J  i  (kr ) 


Pi+  Jo(kr) 


is  an  eigenvector  of  M  if: 


*o  -(a>) ) 


p,  ■ 


to 


ca-i 


p3  =  -iJl  Pi 


c2a2  =  fc2-uM2((£bf-l)c2  . 


i  az 


j  =  1,2, 3,4  (2-12) 


P2=l 


(2-13) 


P4  =  -i^b 


(2-14) 


In  (2-12)  the  index  j  denotes  the  distinct  eigenvectors  corresponding  to 
different  solutions  of  the  dispersion  relation  (2-14)  as  discussed 
below.  In  Equation  13  a  can  be  any  real  number  and  k  can  be  any  positive 
real  number. 

In  (2-14)  we  defined: 

“a  *  f292/4c2 

<*2  -  (Y-l)g2/c2  (2-15) 

.  (u2-(u2 

d  a  b 


By  using  Equation  2-14  we  can  convert  from  the  variable  pair  (k,  a)  to 
(w,  a)  and  we  shall  do  so  later. 


The  eigenvectors  are  normalized  as: 

.  <a,k|k',a'>j  »  N2  i|!±-  6(k-k')  6(a-a')  S.. 

with 


N2 


_ £ _  [.2(a2cV)  +  (a,2-,.2)2] 

(a2cM).2  d'  b 


(2-16) 


(2-17) 


Since  the  dispersion  (2-14)  is  a  quartic  equation  for  <*.,  there  are  four 
orthogonal  eigenvectors,  one  for  each  solution  of  the  quartic: 


U<  . 

1 


(-1)1 

✓2 


[u.2-(k2+a2)c2  +  (-1) 


It1 


((u,2-(k2+a2)c2)2-4k 


2  2 

%c 


2)1/2]V2 


i  «  1,2. 3,4 


(2-18) 


It  can  be  readily  shown  that  for  the  roots  corresponding  to  the  positive 
sign  in  the  square  brackets,  <»..2  _>  u.2  i»l,2  and  for  the  negative  sign 

-J  <  «V  <  -  3-  4. 


D.  THE  INITIAL  CONFIGURATION 

As  the  blast  wave  moves  outwardly  from  a  low  altitude  nuclear 
explosion  It  leaves  behind  most  of  the  energy  of  the  explosion  (some  of 
which  was  initially  radiated  away)  In  the  very  hot  fireball  and  a  region 
surrounding  It  which  is  heated  to  a  much  smaller  degree.  By  the  time  the 
overpressure  reaches  a  value  of  about  .05  the  blast  wave  contains  only 
about  5  percent  of  the  energy  of  the  explosion.  The  general  plan  of  our 
calculation  follow  the  blast  wave  until  the  nonlinear  effects  are 

small  and  jp  „  ve  use  of  linearized  equations;  at  such  point  we  shall 
make  use  or  is  obtained  earlier  which  solve  the  initial  value 
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problem  for  the  linearized  equations.  It  may  appear  that  the  point  at 
which  the  amplitude  of  the  blast  wave  decays  to  .05  would  be  an 
appropriate  time  to  perform  the  coupling  calculation.  Such  is  not  the 
case  however.  As  the  blast  wave  propagates  upwardly  through  the 
atmosphere  it  strengthens  and  lengthens.  A  strengthening  of  the  wave  is 
predicted  by  linear  theory  but  not  very  accurately;  the  lengthening  of  the 
wave  is  an  inherently  nonlinear  effect  which  would  not  occur  at  all  in 
linear  theory. 


A  method  for  calculating  the  nonlinear  propagation  of  the  wave  to 
high  altitudes  was  developed  by  us  and  is  given  in  Reference  8.  The  wave 
assumes  approximately  an  N-wave  shape;  that  is,  it  has  a  shock  wave  at  the 
leading  edge  across  which  the  pressure  jumps  to  some  large  value  followed 
by  an  approximately  linear  decrease  (with  distance  behind  the  shock)  to  a 
value  below  the  ambient  pressure  ending  in  an  abrupt  return  to  ambient. 
Specifically,  we  shall  take  the  hydrodynamic  variables  to  be  given  by 


|f(0)>.  4ir(1 


0 


I 


(2-19) 


where  A  and  L  are  given,  as  a  function  of  the  yield  and  altitude  D  of  the 
event,  in  Reference  8. 


We  have  now,  we  believe,  an  accurate  description  of  the  wave  at 
some  high  altitude  and  It  Is  at  this  point  we  shall  perform  the  coupling 
to  the  freely  propagating  gravity  waves.  At  this  time  the  front  of  the 
wave  Is  no  longer  weak  but  substantial,  possibly  even  considerably  larger 
than  one  for  a  large  low  altitude  explosion.  The  reader  may  wonder  what 
justification  we  have  for  treating  the  problem  In  the  linear  theory  from 
this  point  on.  What  we  offer  is  more  an  excuse  than  a  justification. 
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While  the  amplitude  of  the  wave  is  larger  most  of  the  strength  is 
contained  in  short  wavelengths.  The  content  of  the  wave  at  wavelengths 
which,  as  we  shall  see  later,  are  primarily  responsible  for  coupling  to 
gravity  waves  is  in  fact  small.  Furthermore,  much  of  the  nonlinear 
lengthening  of  the  pulse  has  already  taken  place  (the  reason  for  our  using 
the  method  of  Reference  8  to  propagate  the  wave  to  high  altitudes).  Thus 
the  content  of  the  wave  at  the  most  important  wavelengths  is  small  and 
will  change  only  a  little  for  the  rest  of  the  development.  We  must 
anticipate,  however,  some  dependence  of  our  results  on  the  altitude  at 
which  we  choose  to  do  the  coupling.  Another  reason  for  a  coupling 
altitude  dependence  of  our  results  is  that  the  method  used  in  Reference  8 
calculates  only  an  approximate  wave  form  making  no  attempt  to  calculate 
the  details  of  the  change  in  shape  as  the  wave  propagates.  We  thus 
believe  that  Equation  (2-19)  provides  a  reasonable  initial  condition  for 
linearized  calculation  of  the  amplitude  of  the  gravity  wave  response  and 
obtaining  the  general  features  of  the  wave  form.  Detailed  calculations  of 
the  wave  form  would  require  more  accurate  nonlinear  calculations.  It  is 
perhaps  surprising  that  our  calculated  waveforms  agree  with  the  data  as 
well  as  they  do. 

Since  the  eigenvectors  |k,a>^  form  a  complete  set,  the  vector 
| f (0)>  can  be  expanded  as: 

4  , 

|f(0)>  « l ^  /dk/da  ^  |k,a>i  1<k,a|f(0)>  (2-20) 

To  obtain  the  pulse  at  any  later  time  we  write: 

4  —  ^  (a,  £ 

|f(t)>  ■  /dk/da  ^  e  1  |k,a>j  j<k,a|f(0)>,  (2-21) 

with  <•#■!  given  by  Equation  18. 

For  calculational  purposes  It  Is  advantageous  to  use  (w,  a)  as 
Integration  variables.  The  Integral  (2-21)  can  then  be  written  as 
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I f (t )>  =  {/  du.  J  da 

-U,b  -» 


/  )du,  /  da  } 

r~i  t 

(Dg  /U,  -(jLg 


N2 


ou,a><u.,a  f(0)> 


(2-22) 


where 


Y(*,a) 


2 


(2-23) 


is  the  Jacobian  of  the  transformation.  Note  that  the  sum  over  different 
eigenstates  is  now  represented  by  u.  integrals  over  different  ranges,  the 
first  term  corresponding  to  i=3,4  and  the  third  and  fourth  term  to  i*l  and 
2  respectively. 

The  evaluation  of  «».,a|f(0)>  has  been  given  in  Reference  7.  Here 
we  only  quote  the  final  result 

«t,a|f(0)>  =l^Te'iaD[P*  (—  (sinoO-sino(D-l) )  -  (cosoD+coso(D-L))) 

1  ca2  Lo 

+(ia+kP*)  I  (SI(oD)-SI(a(D-L)))  -  (sinoO+sina(D-L) ) 
o  L 

(cosoO-coso(D-L) ) )]  (2-24) 

ol 


Here 


kz+a‘ 


(2-25) 


and  SI  (x)  Is  the  sine  Integral  function.  In  contrast  to  the  two  layer 
atmosphere  considered  in  Reference  7,  o  is  always  real  in  the  Isothermal 


The  overlap  function  vanishes  as  L+0,  i.e.,  in  the  limit  of  no 
source,  as  it  should.  It  is  also  fairly  easy  to  show  that  <o.,a|f(0)>  is 
finite  for  o-K)  since  the  quantity  in  square  brackets  is  proportional  to  cr2 
in  this  limit. 

From  (2-12)  and  (2-20)  we  can  now  evaluate  the  velocity  in  the 
z-di recti  on, 


Re.  Jdtfc/da  «t.,a| f (0)>e_1  ^t‘az)  Y(u.,a)  (2-26) 

For  the  radial  component  of  the  velocity  we  obtain 

V  =  Bfi.  Jdfc/da  P3(^,a)«*.,a|f(0)>e“i  (*t -az )ili!llA  y(u.,a)  (2-27) 

r  /y  1  N2 


The  corresponding  r,  z-components  of  the  displacement  can  be  found  by 
numerically  integrating  Vz  (t),  Vr  (t)  or  by  closed  form  integration,  e.g. 


d  -  tefdofda  1  <sa|  f  (0)>e“*  (h,t-az)  £oUiO  Y((1(>a) 

Z  VY  u,  1 


/r 


(2-28) 


The  twiddled  quantities  are  related  to  the  observables  through  a 
simple  exponential  factor.  To  make  this  relation  explicit  we  write: 


-ii 

ATT 


(2-29) 


and  we  can  then  show,  assuming  the  scale  height  is  independent  of  z  as  was 
done  in  obtaining  the  linearized  equations.  Equation  (2-1): 


V 

1 


0  exp^)  V 

W  2FT  1 


1  *  z,r 


(2-30) 
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and  —  is  the  relative  overpressure  at  height  D.  Strictly  speaking, 

Equation  2-30  is  incorrect  since  H  depends  strongly  on  z,  varying  by  an 
order  of  magnitude.  Thus  the  absolute  normalization  of  the  pulse  predict¬ 
ed  by  2-30  must  be  modified  to  incorporate,  at  least  approximately,  the 
varying  scale  height.  To  do  this,  we  rewrite  2-30  as: 

V  =  —  D  exp (— Ll2_)  V  ;  i  =  z,  r.  (2-31) 

TP  2Heff  1 

In  our  calculations  we  determine  by  the  following  procedure.  We 
note  that  the  ambient  pressure  at  the  height  z  is 

_  _z _  _  z-D 

P(z)  *  P0e  H(z)  -  P(D)  e  Heff  (2-32) 


Thus: 


^eff 


P(D) 

PTzT 


(2-33) 


where  P(D)  and  P(z)  are  taken  from  tables  of  atmospheric  properties.  For 
example,  for  z«300  km  and  O-lOOkm  we  get  Heff.25km.  A  simpler  estimate 
Is  to  use  Heff  a  1/2  (H(z)  +  H(D))  -  31.4  km,  resulting  in  a  difference 
of  factor  of  2  in  the  overall  normalization. 

To  evaluate  V^,  all  that  Is  needed  to  evaluate  V^.  This  can 
be  done  either  numerically  or  using  any  of  the  available  approximation 

techniques,  such  as  the  stationary  phase.  As  we  shall  see,  only  the  exact 
calculations  are  capable  of  providing  the  detailed  Information  about  the 
time  dependence  of  the  atmospheric  motion,  while  the  stationary  phase 
gives  some  qualitative  features  of  the  TID,  but  does  not  provide  useful 
quantitative  Information. 
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E.  EVALUATION  OF  THE  DOUBLE  INTEGRAL 

The  generic  form  of  the  integral  is  (cf.  Equation  (2-23)) 

V  =  /dou/da  y(a-,a)e_i  ^  (2-34) 

-I  du.  /  da  y(u,,a)e-i  ^  +  (/  d<* +/  du,)  /  da  a  y(u,,a)e  i(i,t  (2-35) 

A  •*  “'a  *  -/J^z 

d 

Since  the  high  frequency  components  of  the  blast  wave  do  not  propagate 
very  far  in  the  atmosphere,  it  appears  reasonable  to  assume  that  the 
second  term  contributes  little  at  large  distances.  We  can  therefore 
approximate  the  integral  by  the  first  term  only.  In  principle,  this  eval 
uation  becomes  very  time  consuming  if  needed  for  a  number  of  different 
values  of  t,z,D,r,L.  To  improve  the  speed  of  the  numerical  calculations 
we  note  that  we  are  interested  in  V  as  a  function  of  time.  Therefore,  we 
construct  tables  of  the  integrals: 

YKT)  -  /da  y(ou,a)e‘Ut  (2-36) 

for  a  fixed  t*T  we  then  have 

V(t)  *  Jd«uY(tu,T)  e1fc(T“t)  .  (2-37) 

Once  Y(^,T)  Is  known  to  any  desired  accuracy,  the  evaluation  of  V(t)  from 
(2-37)  is  quite  rapid.  The  results  of  the  numerical  calculations  are 
reliable  only  for  small  times,  l.e.,  the  first  few  periods  after  the 
arrival  of  the  pulse.  To  calculate  for  later  times  it  is  necessary  to 
develop  approximate  techniques  for  the  evaluation  of  the  double  integral. 

The  method  adopted  Is  that  of  stationary  phase,  which  has  been  used  by 
several  authors  (4).  The  essential  Idea  of  the  stationary  phase  Is  that 
the  major  contribution  to  an  Integral  of  an  oscillatory  function  arises 
from  the  point  at  which  the  phase  is  stationary. 
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To  obtain  a  form  amenable  to  the  stationary  phase  treatment,  we 
note  that: 

Si(x)  =  f (x)sinx+g(x)cosx  (2-38) 

Ju(x)  *  F(x)sinx+G(x)cosx  (2-39) 

where  f,g,  F,G  are  functions  which  we  assume  have  no  oscillatory 
structure.  Strickly  speaking,  this  is  not  correct.  For  example,  in  the 
Bessel  function,  the  argument  in  the  trigonometric  functions  contains 
corrections  of  order  1/x.  However,  since  the  stationary  phase  is  an 
approximation,  this  probably  does  not  cause  large  errors. 

The  phases  in  the  double  integral  are  of  the  form: 

$  =  <*t±kR-a(z-D)±od  ,  (2-40) 

where  d=D  or  d=D-l.  To  obtain  the  stationary  points,  we  need  to  calculate 
the  roots  as,<*.s  of  the  equations. 

U-o  •  W’°  (2-41) 

This  leads  to  8th  order  algebraic  equations  whose  solutions  cannot  be 
given  in  closed  form. 

Since  we  know  that  high  frequencies  do  not  propagate  easily  over 
large  distances,  we  assume  w2<<u,2  .  With  this  assumption  the  dispersion 

(16)  becomes: 

a2  -  -(-4)  +  k2(4)  (2-42) 
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and  the  stationary  points  are 


d.,t?C  ikC 

+  0  .D 


r/c2t2-P2  r 


a  -  ±  *  /t2c2-p2  “  u-  t/r 
s  re  a 


(2-43) 


k  - 


VC 


s  r‘ 


and  a2  *  a2.  Here  C  =  z-D±d  and  p  =  u,  r/u,. 

a  b 


The  stationary  phase  method  then  consists  of  approximating  the  integral 


I  *  /dfc/da  e^^,aV(».,a) 


(2-44) 


by: 


I  -  I 


2* 


S  r a2* 


S.as 


32»(as  ,<■»') 

3d,2 


_  F ( ct.  ,a  )  e'<*<V»s>+"/4> 

V2  s  s 


(2-45) 


and  the  sum  is  over  all  stationary  points. 


After  some  tedious  algebra: 


[ 


3z*(<».,a) 

3a2 


3  ?(as,<».) 


<*<5  ,as  3<*« 


1 


1/2  c  a  i 

■  S  /Oa.  2.2 


W 


(2tV  +a‘r  ) 


1/2 


(2-46) 
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Let : 


Cl  =  z,  C2  3  z-2D,  C3  =  z-L,  c4  =  Z+2D-L 

*!+)s  St'i*  -  ks(5i)R  -  Vi 

-<^s(ci  )t+k s( )R  -  a$c; 


(2-47) 


(2-48) 


and  there  are  8  stationary  phase  points  in  the  integral.  The  evaluation 

of  the  stationary  phase  integral  is  fairly  rapid,  as  compared  to  the  exact 

calculations.  However,  the  stationary  phase  calculation  is  possible  only 

at  time  t>P/c,  and  p>r,  i.e.  after  the  pulse  arrival.  For  y  =  1.4, 

p  =  1.108r  and  the  stationary  phase  derivation  based  on  the  assumption 

o,2<<(».2  is  valid  only  for  t  >  1.5  R/c  by  which  time  most  of  the 
3  ■** 

disturbance  has  died  out.  Nevertheless,  the  stationary  phase  is  useful  in 
making  qualitative  statements  about  the  nature  of  the  oscillatory 
behavior.  For  example,  the  form  of  the  stationary  phase  <t>s  suggests 
that  the  fastest  period  of  the  oscillation  is  about 


as  is  indeed  borne  out  by  our  exact  calculations. 


(2-49) 


Before  proceeding  further  we  discuss  other  stationary  phase 
treatments  found  in  literature  and  point  out  some  drawbacks  of  this 
method.  In  the  exact  calculations  the  velocity  and  displacement  are 
related  by: 

V  -  Re  (/do,  F(u,)eia,t)  (2-50) 


d  -ReUdo^e1*) 


(2-51) 


1 


and  it  is  a  standard  assumption  in  the  literature  that  if  both  V,  d  are 
calculated  from  the  stationary  phase  approximation  to  these  integrals 
then  the.  relation 

d  =  JdtV(t)  (2-52) 

is  satisfied.  We  show  (in  Appendix  A)  that  this  is  not  exactly  correct 
and  it  is  necessary  to  decide  which  integral  is  (better)  approximated  by 
the  stationary  phase.  The  difference  in  evaluating  d  from  (2-51)  or  from 
(2-52)  depends  on  the  structure  of  the  dispersion  relation  and  the  form  of 
the  function  F(u-)  and  the  accuracy  of  the  stationary  phase  approximation. 

A  related  problem  is  the  fact  that  the  stationary  phase  approxi¬ 
mation  does  not  vanish  as  L+0,  i.e.,  in  the  limit  of  no  source.  Both 
these  difficulties  can  be  remedied  by  including  higher  order  terms  in  the 
stationary  phase  expansion,  thus  satisfying  the  relation  to  any  arbitrary 
accuracy.  However,  the  evaluation  of  the  higher  order  terms  involves  some 
extremely  complicated  algebraic  manipulations  (for  the  dispersion  relation 
(2-16)  considered  here)  and  is  somewhat  labor  prohibitive. 

A  standard  stationary  phase  treatment  in  the  literature  is 
carried  out  on  the  Green's  function.5  We  feel  that  such  a  treatment  is 
inconsistent  for  two  reasons.  First,  the  stationary  phase  Green's 
function  is  used  in  an  integral  over  time  which  yields  the  resultant 
displacement  and  velocities.  As  indicated  in  Appendix  A,  this  may  lead  to 
significant,  and  more  importantly  uncontrolled  errors,  thus  diminishing 
the  usefulness  of  such  calculations.  Second,  in  the  Green's  function 
evaluation,  the  stationary  points  are  independent  of  the  source,  contrary 
to  our  results.  A  more  detailed  investigation  is  needed  to  determine  the 
accuracy  of  these  approximate  treatments. 


22 


F. 


NUMERICAL  CALCULATIONS 


To  investigate  the  possible  approaches  to  calculating  the 

atmospheric  gravity  waves  we  have  developed  two  computer  codes.  The  first 

one  evaluates  Vz,  Vr,  dz,  dr  as  functions  of  time  for  fixed  r,  D, 

L  exactly  from  Equations  (2-26)  and  (2-27)  and  includes  also  the  high 

frequency  (u,  >u,  )  contributi ons .  The  second  code  calculates  V  ,  d  ,  V  , 

3  z  z  r 

d^  using  the  stationary  phase  methods.  The  calculations  were  carried  out 
for  a  range  of  radial  distance  r  =  (1000  -  4000)  km  and  yields  Y  ranging 
from  1  kT  to  30  MT  at  observation  height  z  =  300  km.  In  most  calculations 
the  pulse  was  coupled  at  D  =  100  km  and  the  value  of  y  =  1.4  was  used. 

Our  calculations  are  normalized  as  in  Equation  (2-31)  and  c  =  334  m/sec. 


In  order  to  test  our  codes  we  have  carried  out  calculations  in 
the  region  where  the  pulse  originates.  Although  the  calculation  is  quite 
time  consuming  numerically  and  requires  the  inclusion  of  the  high 
frequency  parts  of  the  integral,  we  have  been  able  to  reproduce 
successfully  the  initial  pulse.  We  have  also  tested  the  importance  of 
including  the  high-frequency  contribution  to  the  mode  sum.  For  large 
radii  (r  >  1000  km),  this  contribution  turns  out  to  be  negligible  for  all 
times  t>R/c,  R2  =  rz+zz. 

The  next  calculation  we  wish  to  discuss  is  the  comparison  of  the 
stationary  phase  and  the  exact  numerical  calculations  illustrated  in 
Figure  1.  As  is  readily  evident  from  these  results,  the  stationary  phase 
is  quite  reliable  in  predicting  the  general  time  dependence  of  the  gravity 
waves  even  at  fairly  early  times.  However,  the  stationary  phase  is 
totally  unreliable  In  predicting  the  amplitude  of  the  gravity  wave 
response,  except  at  very  late  times  (after  3  or  4  periods).  This  is 
unsatisfactory  since  by  this  time  the  disturbance  has  essentially  died 
out.  At  early  times,  the  amplitudes  predicted  by  the  stationary  phase  are 
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t-R/c ( 10H  sec) 

Radial  displacement  dr(t)  plotted  vs.  time  t.  Comparison  of 
exact  and  stationary  phase  calculations  at  z  =  300  km, 
r  *  4000  km,  and  Y  *  1  HT. 


several  orders  of  magnitude  greater  than  the  exact  calculations  and  no 
stationary  phase  calculations  can  be  carried  out  for  t  <  u.ar/u.bc. 

This  suggests  that  in  order  to  understand  the  early  time  behavior  of 
gravity  waves  it  is  necessary  to  carry  out  exact  calculations  and  the 
stationary  phase  is  qualitative  at  best. 

With  this  in  mind,  we  set  out  to  investigate  numerically  the 
properties  of  the  gravity  wave  response  as  functions  of  the  yield,  hori¬ 
zontal  distance  and  y  and  compared  with  observations.  In  the  first  set  of 
calculations  we  calculated  the  displacement  dr  for  r=1000,  2000,  and 
4000  km  (i.e.  r-R)  for  the  yield  Y  =  10MT.  The  results  show  an  initial 
outward  displacement  of  the  atmosphere  of  about  20-30  km,  its  magnitude 
essentially  independent  of  the  horizontal  distance,  followed  by  an  inward 
displacement  whose  amplitude  decreases  with  R  as  1/R.  The  period  of 
oscillations  increases  linearly  with  R,  in  agreement  with  the  stationary 
phase  analysis  and  observations.  These  results  are  illustrated  in  Figure 
2  where  we  plot  dr(t),  t  *  *R/c.  Note  that  the  response  is  essentially 
zero  for  t<1,  i.e.,  prior  to  the  pulse  arrival,  as  it  should  be.  An 
important  information  relevant  to  test  detection  is  the  dependence  of  the 
gravity  wave  response  on  the  yield  of  the  weapon  exploded.  This  depen¬ 
dence  does  not  follow  trivially  from  the  formal  expressions  since  differ¬ 
ent  yields  correspond  to  different  A  and  L,  neither  of  which  have  a  simple 
yield  dependence.  We  have  carried  out  exact  numerical  calculations  for 
yields  from  lkT  to  30  MT.  Some  of  these  calculations  are  plotted  in 
Figures  3  and  4  where  we  plot  dr(t)/Y  for  r  *  1000  and  4000  km  and  Y  = 

,1MT  -  30  MT.  It  is  readily  seen  that  in  the  first  period  dr(t)  is  very 
nearly  proportional  to  the  yield.  The  yield  scaled  curves  begin  to  differ 
at  later  times.  The  linear  dependence  on  yield  remains  valid  at  later 
times  for  the  maximum  amplitudes,  however.  The  absolute  displacements 
dr (t )  are  plotted  in  Figure  5,  for  Y*lkT  and  10  kT  at  r*2000  km.  The 
calculated  yield  dependence  of  the  response  cannot  be  readily  deduced  from 


eld  dependence  of  the  radial  displacement  at  z  =  300  km 
=  1000  km. 


t-R/c(103  sec) 

d  dependence  of  the  radial  displacement  at  z  =  300  km 
4000  km. 


t-R/c ( 103  sec) 

d  dependence  of  the  radial  displacement  at  z  =  300  km 


either  the  exact  equations  or  the  stationary  phase  calculations.  Another 
numerical  analysis  has  been  carried  out  to  study  the  y  dependence.  Our 
results  indicate  that  the  period  of  the  gravity  wave  response  is  approxi¬ 
mately  proportional  to  17/7-1  in  agreement  with  the  stationary  phase 
prediction,  and  the  amplitude  is  roughly  independent  of  y.  This  is  illus¬ 
trated  in  Figure  6  for  Y  =  1.67  and  y  *  2. 

For  the  calculated  radial  distances,  the  vertical  displacement 
dz(t)  is  found  negligible  compared  to  dr.  For  observation  carried  out 
at  and  magnetic  dip  angle  and  4,  the  angle  of  propagation  relative  to  the 
magnetic  meridian,  the  vertical  height  variation  is  then 

h (t )  *  d  (t)  sin29cos<t>/2. 
r 

This  gives  variations  in  height  calculated  at  Tonga  (9=-41.5°,  4>=8.1°,r* 
3900  km)  for  the  megaton  range  Housatonic  test  of  order  20-40  km,  in 
agreement  with  reported  observations.1-5  It  should  be  noted  that  for 
bursts  in  the  1  kT  range,  the  calculated  height  variations  are  several 
meters,  well  below  the  sensitivity  of  the  standard  ionosond  observations 
and  probably  well  masked  by  the  noise  due  to  ionosphere  motion  caused  by 
other  (natural)  mechanisms.  Thus  it  is  not  likely  that  the  gravity  wave 
can  be  used  as  a  reliable  guide  in  determining  the  yield  of  the  weapon 
exploded  for  low  yield  (YclOkT). 

In  all  the  calculations  above  we  had  coupled  the  initial  pulse 
at  0*100  km.  This  Is  an  arbitrarily  chosen  altitude  which  may  be  motivat¬ 
ed  by  the  fact  that  the  atmosphere  is  almost  discontinuous  at  this 
height.  The  calculation  of  | f (0 ) >  using  methods  of  Reference  7  should  be 
reasonable  below  this  height  since  we  expect  the  pulse  to  contain  very 
little  of  low  frequency  components.  We  also  carried  out  calculations  for 
D  *  150  and  200  km.  The  calculations  indicate  a  strong  dependence  on  D 
which  Is  evidently  associated  with  the  almost  discontinuous  behavior  of 
the  atmosphere  around  D  *  100  km,  which  changes  rapidly  the  values  of  A 
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and  L  calculated  in  the  nonlinear  model  of  Reference  8  which  are  used  as 
input  in  present  calculations.  It  should  be  noted  that  | f ( 0 ) >  is  not  a 
solution  of  the  linearized  equations  of  motion  so  that  the  strong 
O-dependence  is  not  entirely  unexpected.  In  all  calculations  of  the 
D-dependence,  the  best  agreement  is  in  the  size  and  position  of  the  first 
peak  and  the  choice  of  D=100km  for  our  other  calculations  is  more 
motivated  by  the  agreement  with  observations  then  rigorous  arguments. 
Whether  this  problem  can  be  corrected  by  using  a  two  layer  atmospheric 
model  as  outlined  in  Appendix  B  or  some  other  means  requires  further 
numerical  work. 


SECTION  3 

RISING  FIREBALL  AS  A  SOURCE  OF  GRAVITY  WAVES 


A.  INTRODUCTION 

Calculations  in  the  previous  section  indicate  that  observed 
gravity  waves  from  large  yield,  low  altitude  nuclear  explosions  can  be 
produced  by  the  source  which  provides  a  high  altitude  blast  wave,  namely 
the  initial  rapid  release  of  energy  at  the  burst  location.  While  it  has 
long  been  recognized  that  nuclear  explosions  produce  distant  gravity 
waves,  the  details  of  the  source  mechanism  and  its  relation  to  amplitude 
and  time  history  of  the  observed  gravity  wave  have  not  before  been  firmly 
connected. 


While  it  has  been  accepted  that  the  energy  source  is  important, 
it  has  been  suggested  by  several  authors  that  the  fireball,  which  is  sub¬ 
sequently  produced  by  the  energy  release,  undergoes  motion  which  may  also 
be  a  significant  source  for  gravity  waves.9’10’11  The  fireball  is  formed 
within  seconds  of  the  burst  and  it  consists  of  a  nearly  spherical  hot 
region  of  greatly  reduced  density  which  is  approximately  in  pressure 
equilibrium  with  the  surrounding  cool  air.  However  due  to  the  ambient 
density  stratification  there  is  a  net  buoyant  force  on  the  fireball  which 
causes  an  initial  upward  acceleration  and  a  subsequent  rise  at  roughly  a 
constant  velocity.  The  complex  details  of  the  fireball  rise  dynamics 
Include  entrainment  of  cooler  air,  vortex  formation  and  expansion  due  to 
the  decrease  in  ambient  pressure  with  Increasing  altitudes.  The  fireball 
motion  does  not  cease  immediately  upon  reaching  its  equilibrium  altitude. 
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rather  it  tends  to  overshoot  this  altitude  and  oscillate  about  it  at 
approximately  the  Brunt  period.  At  the  same  time  that  the  fireball  oscil¬ 
lates  around  the  equilibrium  altitude,  it  also  expands  and  contracts  in 
unison  due  to  stratification  of  the  ambient  pressure. 

In  this  section  we  will  model  several  aspects  of  the  fireball 
motion  in  order  to  find  the  significance  of  such  motion  as  sources  for 
gravity  waves  relative  to  the  blast  wave  source.  This  will  be  done  using 
two  different  models  for  fireball  rise  at  a  constant  velocity  along  with 
models  for  fireball  oscillation  about  the  equilibrium  altitude  and  for 
fireball  radial  expansion  and  contraction  near  the  equilibrium  altitude. 

It  will  be  concluded  that  for  the  liner  theory  used,  while  the 
fireball  motion  does  indeed  act  as  a  gravity  wave  source,  this  source 
mechanism  is  not  important  when  compared  to  the  initial  deposition  of 
energy  which  produces  the  blast  wave  so  that  fireball  motion  can  be 
ignored  for  the  purposes  of  gravity  wave  production. 

B.  BACKGROUND 

A  low  altitude  nuclear  burst  quickly  dumps  energy  into  a  local¬ 
ized  region  which  is  approximately  spherical.  A  portion  of  this  energy  is 
radiated  away  in  the  form  of  waves,  linear  and  nonlinear,  in  both  the 
acoustic  and  gravity  wave  frequencies.  In  the  previous  section  the  grav¬ 
ity  wave  from  the  blast  wave  produced  by  the  initial  energy  deposition  was 
calculated.  However  after  the  blast  wave  has  departed  from  the  immediate 
region  of  the  burst  there  remains  a  fireball  which  undergoes  motion  which 
can  serve  as  an  additional  source  for  gravity  waves.  The  fireball  is 
formed  within  seconds  of  a  nuclear  explosion  and  it  consists  of  a  heated 
region  of  a  debris-air  mixture  which  is  approximately  in  pressure  equilib¬ 
rium  with  the  surrounding  cooler  air.  The  fireball  temperature  is  about  1 
eV  and  the  density  is  about  10%  of  ambient.  The  fireball  volume  is  pro¬ 
portional  to  the  yield  so  the  radius  "a"  scales  like 
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(cm) 


(3-1) 


a  =  10s  (Ymt)1/3 

so  a  1  MT  burst  gives  a  radius  of  about  1  km. 

Since  the  underdense  fireball  resides  in  the  stratified  atmos¬ 
phere,  there  is  a  buoyant  force  which  tends  to  cause  the  heated  region  to 
rise.  The  dynamics  of  fireball  rise  include  the  interaction  of  the  buoy¬ 
ant  force  with  the  distortion  of  the  fireball  which  entrains  cooler  air 
and  forms  a  vortex  ring.  The  net  effect  of  this  complex  behaviour  is  a 
gross  rise  of  the  original  fireball  material  at  a  velocity  which  is  nearly 
constant  and  is  approximately 

V  =  Zip'  (3-2) 

where  g  is  the  acceleration  of  gravity  and  S  is  a  measure  of  density 
defect  of  the  fireball  0  =  ( Pg  -  p)/p0.  This  velocity  is  that  appropriate 
to  thermal  rise  which  was  given  by  Scorer.12 

The  flow  outside  a  couple  of  fireball  radi i  is  approximately 
that  seen  around  an  object  rising  at  constant  velocity.  At  smaller  radii 
the  features  we  see  are  much  more  complex.  However,  since  the  size  of 
this  region  is  small  compared  to  the  wavelengths  of  distant  gravity  waves, 
we  expect  that  the  details  of  the  motion  will  not  be  important  in 
dictating  the  general  character  of  that  gravity  wave. 

Fireball  rise  from  a  1  MT  burst  continues  for  a  few  hundred 
seconds  until  it  reaches  such  an  altitude  that  the  density  of  the  expand¬ 
ing  fireball  equals  that  of  the  ambient  air.  This  stabilization  altitude 
('  20  km  for  a  1  MT  burst)  is  dictated  by  the  yield.  The  fireball  does 
not  immediately  settle  at  the  stabilization  altitude  but  tends  to  oscil¬ 
late  about  it  after  an  initial  overshoot  of  the  order  of  a  scale  height. 
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This  oscillation  takes  place  at  about  the  Brunt  frequency  and  may  continue 
for  a  few  cycles.  At  the  same  time  the  fireball  is  bobbing  up  and  down, 
it  also  undergoes  radial  expansion  and  contraction.  The  motion  of  the 
fireball  about  the  stabilization  altitude  will  act  as  a  source  for  gravity 
wave  generation. 

To  date  there  has  not  been  a  convincing  calculation  of  the 
extent  to  which  fireball  rise  is  a  source  of  gravity  wave.  Pierce10  has 
presented  a  formalism  which  allows  gravity  wave  source  terms  to  be 
included  in  the  equations  in  terms  of  moments  of  an  excluded  region  such 
as  the  fireball.  However  no  examples  have  been  worked  out.  Murphy  and 
Kahalas11  have  discussed  some  features  of  fireball  rise  sources  but  again 
no  calculations  have  been  made.  The  work  of  Tolstoy  and  Lau  ,  although 
restricted,  is  the  only  example  of  an  attempt  to  find  the  quantitative 
relation  between  fireball  rise  and  gravity  waves. 

Tolstoy  and  Lau  have  solved  for  the  distant  effects  of  a  force 
source  moving  upward  at  constant  velocity  in  an  incompressible  medium 
which  is  bounded  by  rigid  horizontal  surfaces  at  z=0  and  z=500  km.  The 
solution  is  expressed  as  the  sum  over  the  infinite  set  of  modes  where  the 
coefficients  are  evaluated  in  terms  of  source  characteristics.  They  find 
that  for  a  1  MT  surface  burst  at  a  range  of  10,000  km  and  an  altitude  of 
250  km  vertical  displacements  of  about  5  meters  will  be  generated.  The 
amplitude  is  found  to  scale  like  (yield)5  ^3.  An  amplitude  of  5  meters 
is  quite  small  compared  to  observed  displacements  of  tens  of  kilometers 
for  large  yield  bursts.  Furthermore  the  time  behavior  of  Tolstoy  and 
Lau's  displacement  indicates  that  a  period  of  about  20  minutes  dominates 
the  behavior.  This  behavior  is  difficult  to  understand  since  the  Brunt 
period  is  *  4  min  so  the  dominant  period  would  be  expected  to  be  Brunt 
Period  *  (R/z)  ~  160  minutes.  The  authors  further  calculate  the  gravity 
wave  effects  from  a  mass  injection  source  (referred  to  as  an  explosive 
source)  which  gives  a  gravity  wave  amplitude  which  scales  like 
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(yield)1/3.  It  is  found  that  for  low-altitude  bursts  the  fireball  rise 
mechanism  dominates  for  yields  greater  than  1  MT.  However  the  authors 
point  out  that  an  energy  injection  source  will  produce  a  gravity  wave 
amplitude  which  scales  like  the  yield  but  no  comparison  is  made  with  the 
fireball  rise  mechanism. 

The  force  source  used  by  Tolstoy  and  Lau  was  drawn  from  work  by 
Warren14  who  calculated  the  force  on  objects  moving  at  constant  velocity 
in  a  stratified  incompressible  inviscid  medium.  The  force  on  the  objects 
results  exclusively  from  the  stratification  of  the  medium.  Warren  found 
that  a  sphere  of  radius  "a"  moving  at  a  velocity  V  in  a  medium  of  density 
p0  and  scale  height  H  requires  a  sustaining  force  of 


where  Rd  is  a  function  of  velocity,  through 


For  a  sphere,  Rrf  is  a  function  which  peaks  at  about  1  for  v  »  0.8  and  is 
proportional  to  v2  for  small  v.  This  force  is  referred  to  as  the  wave 
resistance  because  the  energy  expended  through  it  is  radiated  away  as 
gravity  waves.  It  should  be  noted  that  F  peaks  (v  =  0.8)  when  the  time 
for  the  fireball  to  rise  a  diameter  (2a/V)  is  approxi mately  equal  to  the 
Brunt  period  ( /g/H/2ir  for  incompressible  fluids).  Since  the  buoyant  force 
is  at  most 


8 

max 


4 it  p0  a3/3  , 


(3-5) 


the  wave  resistance  is  generally  less  than  the  buoyant  since  B/F  =  4H/a 
and  H  is  about  5  km  at  low  altitudes. 


Elements  of  Warren's  theory  have  been  verified  by  experiments  by 
Mowbray  and  Rarity. 15  In  these  experiments,  a  sphere  in  a  stratified 
medium  was  moved  at  a  constant  velocity  while  looking  for  evidence  of 
gravity  waves.  It  was  found  that  no  observable  wave  pattern  was  set  up  if 
the  velocity  was  too  great  (in  Warren's  terms  v  <  0.2  gave  no  pattern)  but 
that  a  smaller  velocity  (corresponding  to  v  «  0.8)  gives  rise  to  a  good 
wave  pattern  with  periods  suitably  related  to  the  Brunt  period. 

Other  experiments  by  McLaren,  Pierce,  Fohl  and  Murphy16  are 
relevant  to  the  generation  of  gravity  waves.  These  experiments  consisted 
of  the  release  of  underdense  fluid  elements  in  a  stratified  medium  with 
subsequent  observations  of  fluid  motion  well  away  from  the  source.  The 
results  indicate  no  gravity  waves  are  generated  during  the  rise  phase  of 
the  motion.  However  gravity  waves  are  generated  following  the  rise  when 
the  sample  fluid  element  oscillates  approximately  two  times  about  its 
stabilization  height.  It  is  not  clear  that  the  failure  to  generate  waves 
during  the  rise  is  in  contradiction  to  Warren's  theory  since  inadequate 
information  is  provided  to  find  a  value  of  v  for  comparison.  Furthermore 
it  is  not  obvious  that  a  solid  body  under  an  external  force  should  gener¬ 
ate  the  same  waves  as  a  miscible  fluid  element  moving  at  the  same  velocity 
under  strictly  hydrodynamic  forces.  In  any  case,  the  results  of  McLaren 
et  al.,  indicate  that  the  oscillations  about  the  stabilization  altitude 
should  be  considered  as  a  gravity  wave  source  as  an  element  of  the  fire¬ 
ball  rise  mechanism. 

C.  GRAVITY  HAVE  THEORY 

The  hydrodynamic  equations  of  motion  of  a  perfect  gas  in  the 
presence  of  sources  can  be  written  as 
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Where  0  is  a  mass  source  (sec-1),  F  is  an  acceleration  source  (cm/sec2) 
and  S  is  a  heat  source  (ergs/gm/sec) .  When  these  equations  are  linearized 
in  an  exponential  atmosphere  of  scale  height  H  with  constant  sound  speed 
c,  the  Fourier  transform  of  the  relative  pressure  perturbation  can  be 
written  as 
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=  M 


where  w.  is  the  Brunt  frequency  and  w  is  the  acoustic  cutoff  frequency 
d  a 

defined  by  Equation  (2-15). 
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The  formal  solution  for  the  pressure  pulse  resulting  from  a  source  M  can 
be  written  as 


=  eZ/2H  /  dr1  e'Z'/2H  /  di  G(ir-r'|,  t)  M(r',  w)  (3-8) 
P  1  1 

where  G  is  the  Green's  function  for  the  pressure  equation.  In  general 
there  is  no  closed  form  solution  for  the  time  domain  Green's  function 
However  the  behavior  in  the  limit  of  low  frequencies,  u.  «  is  readily 
expressed  as 


G  |p  =  -  J0  (*  (t2-(R/c)2) l/2)  9(t-R/c)  (3-9) 

Lr  4irR  C 

where  the  characteristic  late  time  frequency  is 


for  a  source  to  observation  range  of  R  and  sound  speed  c.  The  function  0 
is  a  unit  step  function  which  allows  no  response  until  t  >  R/c.  While  the 
low  frequency  Green's  function  will  not  accurately  reproduce  behavior  for 
times  near  R/c,  at  observation  locations  for  which  z  «  r  the  later  time 
behavior  is  well  represented.  This  has  been  verified  in  a  previous 
report.6  As  a  result,  since  we  are  interested  in  behavior  at  points  at 
thousands  of  kilometer  ground  range  from  the  source,  but  at  ionospheric 
heights,  we  feel  that  the  low  frequency  limit  is  a  useful  tool  for  explor¬ 
ing  the  implications  of  different  source  models. 


The  source  function  M,  in  the  low  frequency  regime  can  be 

♦ 

expressed  in  the  time  domain  as  a  rather  simple  function  of  0,  F  and  S 
when  it  is  noted  that  i  u,+-  Thus  for  <*,  «  we  find 
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A  useful  example  of  a  low  frequency  response  is  provided  by  a 
source  which  provides,  in  the  acoustic  range,  an  N-wave  pressure  fluctua¬ 
tion.  As  shown  in  that  previous  report17  such  a  source  of  temporal  length 
To  gives  a  pulse  of  the  form 
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where  N(t)  is  a  unit  amplitude  N-wave  and  the  normalization  A  can  be  fixed 
in  terms  of  the  N-wave  amplitude  at  any  location.  Thus  if  the  N-wave  has 
an  amplitude  (5p/p)o  at  a  vertical  location  z0,  the  low  frequency  gravity 
wave  associated  with  the  same  source  will  have  a  peak  value  of 
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for  a  time  t  just  after  R/c.  Note  that  (i£.)oin  this  equation  is  the  same 
as  in  Equation  (2-30).  This  example  will  serve  as  a  useful  standard 
against  which  to  compare  the  effects  of  other  sources. 
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D.  RISING  FIREBALL  MODELS 

In  this  subsection  we  will  present,  in  the  context  of  gravity 
wave  theory,  solutions  for  a  set  of  source  models  which  simulate  several 
aspects  of  fireball  motion.  These  gravity  waves  will  be  compared  with 
those  from  a  blast  wave  source  in  order  to  establish  the  dominant  source 
mechanism  for  gravity  waves  from  nuclear  explosions.  The  models  include 
two  which  provide  fireball  rise,  one  which  simulates  the  vertical  bobbing 
motion  about  the  stabilization  altitude  and  one  which  describes  the  coher¬ 
ent  radial  oscillations,  or  breathing,  while  the  fireball  oscillates 
around  the  stabilization  altitude. 

D.l.  Moving  Source  of  Force 

Initially  we  will  consider  the  source  model  used  by  Tolstoy  and 
Lau.  This  is  a  point  source  of  force  in  the  vertical  direction  which 
moves  at  a  constant  velocity  V  for  an  interval  of  time  T.  The  accelera¬ 
tion  per  unit  volume  imposed  is 

F  =  —  *  «(x)  «(y)  «(z-Vt)  [9ft)  -  Oft-T) ]  (3-14) 

2  P 

where  0  is  the  unit  step  function.  HeretF'will  be  taken  as  the  net  force 
on  a  sphere  moving  at  uniform  velocity  through  a  stratified  medium  as 
given  by  Warren.  The  source  can  be  considered  as  a  point  because  the 
gravity  wave  wavelengths  of  interest  greatly  exceed  the  fireball  size. 

The  low  frequency  source  mlf  contains  a  term  for  a  force  source 
which  is  proportional  to  the  second  time  derivative  as  well  as  the 
z-derivatlve  of  this  expression.  Since  the  source  region  is  small  and  the 
force  has  a  single  sign,  we  can  neglect  the  z-derivative  contribution. 
Therefore  the  low-frequency  pressure  fluctuations  from  this  source  mechan¬ 
ism  is 
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The  three  terms  lead  to  contributions  which  are  of  the  form 
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where  R  =  [r2  +  (z-z')2]1/2.  Note  that  the  third  term  begins  at  zero  at 
t  =  R/c  and  grows  linearly  in  time  at  early  time  while  the  other  two  terms 
have  their  peak  (although  this  is  not  completely  obvious)  at  t  immediately 
following  R/c.  The  first  two  terms  are  obtained  by  evaluation  of  the 
difference  of  the  expressions  at  the  indicated  values  of  z'  and  t'.  If 
t<  R/c  +  T,  then  only  the  upper  values  are  used  since  the  second  step 
function  to  cut  off  the  source  should  not  be  invoked. 


Now  an  upper  bound  on  the  contributions  of  the  three  terms  can 

be  found  by  evaluating  the  first  two  at  t  =  R/c  and  the  third  at  t=R/c  +T 

when  T  *  ou-1.  This  gives 
c  3 


For  the  parameters  considered  later  for  distant  observations  of  a  1  MT  sea 
level  burst,  it  will  be  found  that  the  second  term  is  the  dominant  element 
of  the  expression. 

D.2  Rising  Mass  Dipole  Source 

The  flow  of  air  around  a  rising  fireball,  at  points  outside  the 
fireball  itself,  resemble  a  dipole  flow  around  a  moving  source  pair 
(source  and  sink)  separated  by  roughly  a  fireball  diameter.  This  combina¬ 
tion  of  motion  and  point  sources  can  easily  be  simulated  to  provide  an 
alternate  model  for  the  fireball  rise  mechanism  as  a  gravity  wave  source. 
Since  it  is  not  clear  exactly  which  aspect  of  the  fireball  rise  is  most 
important  in  production  of  gravity  waves,  it  is  best  to  provide  models  for 
several  features. 
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A  pair  of  rising  mass  sources  displaced  vertically  by  a  separation 
A z  and  acting  over  a  finite  time,  T,  can  be  described  by 


0  S(x)  5(y)  [fi(z-Vt)  -  S(z-Vt-Az)]  [  9  (t )  -  9ft-T)  ]  (3-18) 

p  dt 

where  the  net  emission  of  mass  is  exactly  zero.  A  single  one  of  the  pair 
of  sources  emits  mass  at  a  rate  dM/dt  while  the  other  absorbs  mass  at  the 
same  rate.  Since  the  net  effect  of  the  motion  of  the  fireball  (at  con¬ 
stant  velocity  V)  is  to  move  air  from  above  to  below  the  fireball  the 
displaced  mass  will  give 

—  =  pira2V  (3-19) 

dt 

where  "a"  is  the  fireball  radius.  We  will  take  dM/dt  to  be  constant  and 
simply  treat  the  density,  p,  as  a  parameter. 

We  will  now  sketch  the  calculation  required  for  a  single  source 
and  take  the  sum  of  two  following  the  final  step.  Consider  the  source 

Q  -  -  ^  «(x)  «(y)  S(z-Vt-Zo)  [9(t)  -  9(t-T)]  (3-20) 

p  dt 


In  the  low  frequency  limit  the  pressure  response  to  this  source  is 
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After  the  trivial  integrations  the  result  is 
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If  t  <  R/c+T  the  first  term  is  evaluated  only  at  z'=zo  and  the  upper  limit 
of  the  integral  is  z0+V(t-R/c).  The  peak  value  of  the  above  expression 
can  be  estimated  by  evaluation  of  the  first  term  for  t  -  R/c  and  for  the 
second  at  t  *  R/c  +T.  This  gives  (for  <*c  «  u<b  and  u^T  «  1) 
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The  corresponding  result  for  two  sources  of  opposite  sign  separated  by  Az 
is  then 
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D.3  Oscillating  Sphere 

The  experimental  work  by  McLaren  et  al  suggests  that  a  buoyant 
object  oscillating  about  its  equilibrium  altitude  acts  as  a  source  of 
gravity  waves.  In  order  to  model  the  contribution  of  fireball  oscilla¬ 
tions  to  gravity  wave  generation  we  consider  a  force  source  which  oscil¬ 
lates  in  sign  and  location  in  such  a  manner  that  the  associated  accelera¬ 
tion  is  approximately  that  of  a  fireball.  A  damped  oscillation  can  be 
represented  by 


Fz  =  A  sin  (o,Qt)  e~at  S(x)  6(y)  5(z-z0-Az  sin  o^t )  9(t)  (3-25) 

where  u,  is  the  oscillation  frequency  and  Az  is  the  amplitude.  If  a 

sphere  of  radius  "a"  oscillates  at  a  frequency  o, ,  the  peak  acceleration 

will  be  a,2  Az  over  a  volume  of  4ira3/3.  Thus  at  a  density  p,  the  net  force 
o 

required  to  impose  such  an  oscillation  will  be 

if  =  u2  Az  4ira3p/3  =  Ap  (3-26) 

o 

The  pressure  fluctuations,  in  the  low  frequency  limit,  can  be 

written  as 
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Now  for  distant  observation  points  u,  is  a  small  compared  to  a  since  the 
fireball  motion  is  damped  out  after  a  few  oscillations  at  approxi mately 
the  Brunt  frequency.  That  is 


hj  «  a  * 
c 


Vn 


(3-28) 


where  n  is  the  number  of  oscillations  before  damping.  As  a  result  the 
bulk  of  the  time  variation  of  the  integral  is  provided  by  the  exponential 
factor.  This  leads  to 
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Thus  for  the  above  restriction  on  a,  the  peak  fluctuation  will  be 


z-zq 

2h 

.  e  w.  (J.3  Az  4ira3 

iE.)  - _ *L_i! _  (3-31) 

p  max  12wRg2H( y-l ) 

D.4  Breathing  Sphere 

As  the  fireball  bobs  vertically  about  its  equilibrium  position, 
it  also  undergoes  expansion  and  contraction  cue  to  the  associated  varia¬ 
tion  of  ambient  pressure.  While  there  is  no  direct  evidence  to  suggest 
this  is  an  important  gravity  wave  source,  it  may  be  that  this  mechanism 
produced  the  waves  seen  experimentally  by  McLaren  et  al  since  they  could 
not  separate  the  effects  of  bobbing  and  breathing.  Therefore  we  will  now 
give  a  model  for  this  breathing  which  is  appropriate  as  a  gravity  wave 
source. 
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The  breathing  of  the  fireball  can  be  simulated  by  an  oscillating 
mass  injection  which  provides  no  net  mass  change  over  a  cycle.  We  will 
take  the  mass  change  to  be 


AM  ( t )  =  AM  sin  ( u,  t ) 
o  o 


(3-32) 


It  operates  for  a  single  cycle,  is  continuous  and  acts  at  a  point.  The 
source  for  our  model  is  then 


0  =  —  am  6(r)  cos  (u,  t) 


0  <  t  <  2m/u. 


(3-33) 


so  the  low-frequency  response  is 
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The  terms  in  the  square  bracket  have  a  time  scale  -  <»<■*  while  Jo  has  a 
time  scale  -  Consequently,  since  o.  ■  J0  can  be  expanded 

about  t  -  t  and  the  integral  can  be  directly  evaluated.  When  this  is 
done,  one  finds  that  the  first  three  terms  in  the  series  give  no 
contribution  to  the  integral  so  that  the  result  is 
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so  the  peak  value  at  t  =  R/c  will  be 


<Spj 
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A3) 


E.  EFFECTS  FROM  FIREBALL  MOTION  MODELS 


(3-361 


In  order  to  find  the  magnitude  of  the  gravity  wave  response  to 
the  fireball  motion  models  presented  in  subsection  0,  it  is  convenient  to 
use  a  consistent  set  of  parameters  appropriate  to  a  typical  case  of  inter¬ 
est.  In  this  we  will  follow  Tolstoy  and  Lau  who  considered  the  example  of 
a  one  megaton  low-altitude  explosion  as  observed  at  a  distant  point  in  the 
i onosphere. 

We  will  use  an  observation  height  of  z  =  250  km  where  the  dens¬ 
ity  is  p  -  6  x  10- 14  gm/cm3.  When  this  is  compared  to  the  sea  level  value 
( p0  =  1.2  x  10"3),  it  is  seen  that  a  mean  scale  height  of  H  »  10.5  km  is 
required  to  provide  the  proper  ratio  for  an  isothermal  atmosphere.  Fire¬ 
balls  are  initially  at  about  10%  ambient  density  so  the  density  defect  is 
8  -  .9.  For  a  1  MT  burst  the  initial  fireball  radius  is  a  =>  1  km.  There¬ 
fore  the  fireball  rise  velocity,  for  g  =  950  cm/s2,  is 


V  *  /agB  ■  9  x  io3  cm/s  (3-37) 

The  force  which  must  be  exerted  on  a  sphere  to  sustain  this  motion  in  a 
stratified  medium  is  then 


This  is  an  upper  limit  since  the  peak  value  as  a  function  of  velocity  has 
been  taken.  (Note  that  this  net  force  is  much  less  than  the  buoyant  force 
4/3na  p<j.)  Tolstoy  and  Lau  indicate  that  the  rise  time  of  the  fireball 
before  reading  its  equilibrium  altitude  is 


T  »  300  sec 


This  is  consistent  with  experiment  by  McLaren  et  al  which  suggest  the  rise 
takes  about  one  Brunt  period.  In  this  example 


CL 


b 


1.6  x  10-2  sec-1 


so  the  Brunt  period  is  =  2*/^  »  390  seconds.  To  complete  the  para¬ 
meters  required  we  use  Y  =  1.4,  observation  radius  of  r  =  104km  *  R  giving 
cl^  =  4  x  io-4  sec"1  and  c  =  374  m/s. 

Now  the  rising  force  model  of  Equation  (3-18)  gives 

— )  *  3  x  10"4  _  Rising  Force 

p  max 

The  result  for  the  rising  dipole  mass  source  from  Equation  (3-25)  for  a 
separation  Az  =  2a  is 

=  4  x  io-3  _  Rising  Dipole  Mass 

p  max 

The  oscillating  sphere  source  provides  the  result  of  Equation  (3-31)  which 
for  an  amplitude  Az  of  5  km  and  a  frequency  of  clq  =  gives 

*  4  x  io-3 

p  max 


Oscillating  Sphere 


Finally  the  breathing  sphere  yields  for  AMq  =  0.1  •  4/3wa3p,  correspondi ng 
to  a  10%  volume  fluctuation. 


£P)  =  2.5  x  10-5  _  Breathing  Sphere 

p  max 

For  comparison  with  the  blast  wave  source  as  provided  in  Equation  (3-14) 
we  take  the  1  MT  values  initialized  at  z0  =  100  km;  (5p/p)0  =  2.1  and 
To  =  100  seconds.  This  gives 

— )  =  1.5  -  Blast  Wave 

p  max 

For  comparison  with  observations  it  is  often  useful  to  express 
fluctuations  in  terms  of  vertical  displacement.  Generally  there  is  not 
simple  relation  between  6p/p  and  displacement  but  in  the  case  of  low  fre¬ 
quencies  a  useful  result  can  be  given  because  in  this  regime  the  motion  is 
nearly  incompressible.  Therefore  pressure  perturbations  result  from 
vertical  displacements  of  the  neutral  air  such  that 

^  -  if.  (3-39) 

P  H 

However  observations  in  the  ionosphere  are  of  vertical  displacement  of  the 
plasma  which  tends  to  move  with  the  component  of  the  neutrals  along  the 
earth's  magnetic  field.  Furthermore  the  gravity  wave  particle  motion 
tends  to  be  radial  from  the  source  so  that  combining  this  information  we 
find  a  neutral  radial  displacement  of 

6r  -  5z  -  -  (— )H  li  (3-40) 

z  z  p 
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giving  a  vertical  plasma  displacement  of 

Az  s  <5r  sin  9^  cos  9d  cos  (3-41) 

where  9^  is  the  magnetic  dip  angle  and  4>  is  the  angle  of  propagation  rela¬ 
tive  to  the  magnetic  meridian.  Thus  for  typical  magnetic  angles  the 
product  of  the  trigonometric  functions  is  '  1/4.  The  vertical  displace¬ 
ment  of  the  ionospheric  plasma  from  a  low  frequency  pressure  fluctuation 
is  then 


.  Sp  .  rH 
p  4z 


(3-42) 


The  blast  wave  calculation  for  1  MT  as  observed  at  10,000  km,  from  previ¬ 
ously  stated  parameters,  will  be  ~  150  km  which  is  about  an  order  of 
magnitude  greater  than  that  found  in  the  previous  section.  This  differ¬ 
ence  is  largely  due  to  different  choices  of  effective  scale  heights  (10.5 
km  versus  25  km).  Use  of  25  km  leads  to  a  reduction  of  about  a  factor  of 
30  or  to  a  net  Azof  5  km  in  the  current  calculation.  This  is  comparable 
with  the  result  of  Section  2.  This  illustrates  a  difficulty  of  making  a 
direct  comparison  with  experimental  results  using  an  isothermal  atmosphere 
model . 


F.  SUMMARY 

The  calculations  of  Section  2  indicate  that  the  gravity  waves 
with  periods  of  up  to  hours  observed  at  thousands  of  kilometers  from  low- 
altitude  large  yield  nuclear  explosions  are  generated  by  the  blast  wave 
mechanism.  In  this  section  we  have  explored  the  possibility  that  the 
fireball  rise  and  oscillation  may  also  make  a  significant  contribution  to 
gravity  waves.  To  this  end  a  set  of  four  models  of  fireball  motion  have 
been  constructed  and  there  effects  calculated  using  the  low-frequency 
theory  in  conjunction  with  point  sources.  It  is  found  that  each  of  the 
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source  models  generates  only  very  weak  gravity  waves  relative  to  the  blast 
wave  mechanism  giving  pressure  fluctuations  which  are  about  three  orders 
of  magnitude  less  than  that  from  the  blast  wave  source.  The  models  used 
for  fireball  motion  are  rather  simple.  However  in  each  instance  approxi¬ 
mations  used  were  chosen  to  maximize  the  effects  generated.  As  a  result 
we  can  be  fairly  confident  that  use  of  more  realistic  models  will  not 
alter  the  basic  result  and  such  calculations  are  not  justified  since  it 
seems  clear  that  fireball  rise  is  not  a  significant  source  of  gravity 
waves  from  low  altitude  nuclear  explosions. 


SECTION  4 

SUMMARY  AND  CONCLUSIONS 


In  the  present  contract  DARPA/AFOSR  tasked  MRC  to  study  the  gene¬ 
ration  of  ionospheric  gravity  waves  by  nuclear  surface  bursts  to  determine 
what  aspects  of  the  explosion  are  the  most  important  source  of  ionospheric 
gravity  waves:  special  consideration  has  been  given  to  the  question  of 
whether  or  not  the  traveling  ionospheric  disturbance  observed  by  Lewis 
Duncan  at  Arecibo  on  22  September  1979  could  have  been  caused  by  a  small 
low  altitude  explosion  and  to  the  more  general  question  of  whether  or  not 
ionospheric  gravity  waves  are  likely  to  be  useful  for  detecting  small, 
low-altitude  nuclear  explosions  at  great  distances  from  the  burst.  We 
have  made  the  necessary  calculations  and  have  obtained  good  agreement  with 
the  data  -  the  first  time  such  agreement  has  been  obtained  from  a  theore¬ 
tical  calculation.  The  calculated  gravity  wave  signal  agrees  with  obser¬ 
vation  in  magnitude,  scales  correctly  with  observation  distanced .e. ,  its 
period  increases  linearly  with  distance,  it  maximum  amplitude  decreases  as 
1/R)  and  with  yield  (i.e.,  its  period  independent  of  yield  and  the  ampli¬ 
tude  proportional  to  yield). 

In  the  semiannual  report6  we  presented  the  results  of  the 
generation  of  ionospheric  gravity  waves  by  the  high  altitude  hot  spot 
above  a  low  altitude  nuclear  explosion;  that  work  primarily  considered 
waves  traveling  In  the  ducted  modes.  In  the  present  report  we  have  calcu¬ 
lated  the  generation  of  gravity  waves  by  the  blast  wave  and  by  the  rising 
fireball.  In  this  work  we  have  considered  isothermal  atmospheres  and  so 
have  only  freely  propagating  modes  (the  coupling  of  the  blast  wave  to  the 
ducted  gravity  modes  in  a  two  layer  atmosphere  was  calculated  in  Reference 
7). 


On  the  basis  of  our  research  we  have  reached  the  following 
conclusions:  1)  the  major  source  of  the  large  amplitude  ionospheric 
gravity  wave  observed  following  the  U.S.  and  Soviet  tests  (the  wave  with 
periods  ranging  from  30  minutes  to  more  than  two  hours  is  the  blast  wave 
and  the  energy  is  transmitted  by  the  freely  propagating  modes;  2)  the 
source  of  the  15  minute  period  wave  that  were  observed  to  go  all  the  way 
around  the  world  is  the  high  altitude  hot  spot  and  the  energy  is  traveling 
in  the  ducted  modes;  3)  the  motions  associated  with  the  rising  fireball 
give  a  smaller  signal  in  the  freely  propagating  modes  than  the  blast  wave; 
4)  it  appears  that  ionospheric  observations  will  be  of  limited  use  in 
detecting  small  (>10kT)  low  altitude  detonations  at  large  distances 
(>1000km)  from  the  burst. 
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APPENDIX  A 

TIME  INTEGRATION  IN  THE  STATIONARY  PHASE 


In  this  appendix  we  illustrate  the  failure  of  the  stationary 
phase  in  reproducing  exactly  the  relation 

d  =  Jdt  V(t)  (A-l) 

from  the  stationary  phase  evaluation  of  the  integrals 

V  .  /d*F(»)  e’(“-t-kr>  (A -2) 

d  »  (A-3) 

To  illustrate  the  source  of  the  problem  we  assume  a  simple  dispersion 
relation: 

u.2  *  m2+k2c2  (A— 4 ) 

Then  the  stationary  point  is  obtained  by  solving 


11  «  0 


3(b 

which  gives: 


*  ct 


(A— 5) 


<b 

s 


mtc 

/tV7 


,  tc>  r 


(A -6 ) 
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and  no  solution  exists  for  tc<  r. 


and 


The  stationary  phase  evalution  of  the  integral  requires 

ill  =  m2r  .  (c2t2-r2)3/2 

3u,2  a.  »  dj  c(n-2-m2)3/2  mr2c 


( A  -  7 ) 


mr/c 

/t2c2-r2 


( A -8 ) 


<j>s  =  m(t2-r2/c2)1/2/c 

The  stationary  phase  approximation  then  gives: 

V  =  F(„  ,  e'*s 

s  2mr2ir  s 

d  .  [(tV-r2)3/2]-1/2  fK)  eHs 

s  2mr2*  ids 


( A -9 ) 


(A-10) 


(A— 1 1 ) 


To  carry  out  the  comparison  further,  we  now  choose  the  form  of  the 
function  F(d).  In  particular,  we  take 


Thus 


F(d) 


[(t2c2-r2)3/2]1/2  c 
2mr2* 


i^c 


ce 


» 


(A-12) 


( A- 13) 
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( A- 14 ) 
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^V-r2 

i  mt  e 


To  check  we  evaluate 


ce  [l  + 


i  4>sc  t 


1  *  V 

Z.  2J  S 


(A-15) 


The  second  term  in  (A-15)  is  different  from  zero  for  all  finite  times 
t >r/c  and  indicates  the  magnitude  of  the  error.  To  improve  the  agreement 
between  /Vsdt  and  ds  it  is  necessary  to  include  higher  order 
corrections  to  the  stationary  phase.  Though  possible,  this  becomes 
quickly  impractical,  in  particular  for  a  complicated  dispersion  relation 
of  the  type  needed  to  analyze  the  gravity  waves. 

The  result  of  this  appendix  indicates  that  the  use  of  stationary 
phase  in  the  analysis  of  the  gravity  waves  is  somewhat  questionable  if  a 
detailed  behavior  of  the  atmospheric  motion  is  to  be  derived.  Our  main 
hope  is  that,  at  best,  the  stationary  phase  is  capable  of  providing 
qualitative  information  about  general  trends,  maximum  displacement  and 
other  gross  features  of  the  physics  following  an  atmospheric  burst.  It  is 
interesting  to  note  that  at  t2>>r2/c2. 


V  ce 


i  mt 


c 

im 


imt 

e 


and  the  stationary  phase  is  exact.  It  therefore  seems  reasonable  to 
conjecture  that  the  stationary  phase  becomes  exact  as  w^t^kr  since  “-s 
approaches  a  constant  at  large  t  and  it  is  possible  to  satisfy  this 
condition. 
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APPENDIX  B 

GRAVITY  WAVES  IN  TWO-LAYER  ATMOSPHERE 


We  now  wish  to  consider  the  gravity  wave  response  propagating 
through  the  freely  propagating  modes  in  a  two  layer  atmosphere.  We 
choose,  arbitrarily,  z=0  as  the  altitude  of  the  discontinuity.  In  each 
layer,  the  eigenvectors  |sa>  given  in  Equation  (2-12)  are  solutions.  We 
denote  by  |u.,a>i  the  solution  in  the  lower  layer  and  by  |u',a>2  the 
solution  in  the  upper  layer.  For  a  disturbance  originating  below  z=0, 
there  will  be  an  upwardly  propagating  wave  ju.,a>i  and  a  reflected  wave 
|<S-a>!  in  the  lower  layer.  In  the  upper  layer,  only  the  outwardly 
propagating  wave  occurs.  Thus  the  general  solution  in  a  two  layer 
atmosphere  is: 


In  terms  of  the  coefficients  defined  in  Section  2,  the  matching  conditions 
at  z=0 

cia2|z=0-  =  c2a2|z=0+  (B'2) 


(B-3) 


(B-4) 
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(B-5) 


AP  i  (a i )+BP i (-a i )+  —  P4(ai)  =  C(P1(a2)+  —  P4{a2)) 

Y-l  Y-l 

Here  a-j  is  used  to  denote  the  wave  propagating  below  (i=l)  and  above 
(i=2)  z=0,  Pj(o,,a)  are  defined  in  Section  2B. 

The  matching  conditions  are  readily  solved  for  A/C,  B/C  and  we 

find: 

A/C  =  [P  1(a2)-P  !  ( -a  ! )  +  JL  (P4(a2)  -  P4(a1))]/(Pl(al)-Pl(-a1))  (B-6) 

Y-l 

B/C  =  [P1(a2)-P1(a1)+  —  (P4 (a 2 )  -  P4 (a x ) )  ]/(P L (a L  )-P l (-a  l ) )  (B-7) 

Y-l 

and  the  normalization  requires 

| A | 2  +  |B | 2  +  |C | 2  =  1  (B-8) 

These  coefficients  together  with  (B-l)  provide  a  complete  set  of  states 
and  we  can  write: 

|f(t)>  =  /d<*-/da |  <sa>  <u-,a|f(0)>e  (B-9) 

and  we  can  carry  out  calculations  as  in  Section  2. 
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